Equilibrium and nonequilibrium applications of lattice-gas models in electrochemistry 
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We discuss applications of statistical-mechanical lattice-gas models to study static and dynamic 
aspects of electrochemical adsorption. The strategy developed to describe specific systems includes 
microscopic model formulation, calculation of zero-temperature phase diagrams, numerical simula- 
tion of thermodynamic and structural quantities at nonzero temperatures, and estimation of effec- 
tive, lateral interactions. We briefly review earlier work, including studies by Monte Carlo simulation 
of the adsorption of urea on Pt(100). As an illustrative example, we discuss in some detail recent 
applications to the underpotential deposition of Cu with sulfate on Au(lll). Experimental and 
numerical results are included for equilibrium coverages and voltammetric currents at slow potential 
sweep rates, and new results from dynamic Monte Carlo simulations are presented. In particular, 
we present simulated results for current transients far from equilibrium, following sudden changes 
in the electrode potential. 
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I. INTRODUCTION 



The recent confluence of electrochemistry and surface science has occasioned a two-way exchange of experimental 
and theoretical techniques between these previously disparate disciplines. Here we discuss electrochemical applications 
of a theoretical method adapted from surface science: computational, statistical-mechanical lattice-gas modeling. In 
Sees. 2 and 3 we briefly review some earlier equilibrium and near-equilibrium work including an application to 

the electrosorption of urea on Pt(100) §-§, and we briefly summarize a strategy to obtain effective, lateral adsorbate- 
adsorbate interactions by comparing model results with experiments performed at or near equilibrium. More extensive 
reviews of this work can be found in Refs. 

As an illustration of the approach, in Sec. 4 we give a "case study" of a particular system: the underpotential 
deposition (UPD) of Cu with sulfate on Au(lll) P,pd| p"8| . First we review the results of a recent Monte Carlo (MC) 
study of equilibrium coverages and voltammetric currents at very low potential sweep rates |]l7|| . This discussion is 
followed by an extension of the model to study its dynamical behavior following a sudden change in the electrode 
potential. In particular, we present new simulation results for the time-dependent coverages, voltammetric currents, 
and microscopic adsorbate structures. The reasonable agreement with recent current-transient measurements by 
Holzle et al. |19) is encouraging. 



II. LATTICE-GAS MODELS OF SUBMONOLAYER CHEMISORPTION 



The models we discuss are defined by the two-component lattice-gas Hamiltonian (l]-||,@-||,|l7j 1 
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Here c^g{0,l} is the local occupation variable for species X; Y^Uj) an d run ov er all nth-neighbor bonds and over 
all adsorption sites, respectively; denotes the effective lateral XY pair interaction through an rtth- neighbor bond; 
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J2 n runs over the interaction ranges; and H.3 contains multi-particle interactions. The change in electrochemical 
potential when one X particle is adsorbed is —fix- Particular systems differ in their binding-site geometries and the 
values and ranges of the lateral interactions. The thermodynamic density conjugate to fix is the surface coverage, 

e x = ^- 1 £ c ^ (2) 

i 

where N is the number of surface unit cells. The mapping between the electrochemical potentials in the lattice-gas 
Hamiltonian, fix, and the bulk activities [X] and electrode potential E is accomplished through the relation, 

fix = Mx + RThx ([X]/[X] ) - zxFE , (3) 

where R is the molar gas constant, T is the absolute temperature, F is Faraday's constant, zx is the effective 
electrovalence of X, and /i^and [X]° are reference values which include the local binding energies. 

The interactions in Eq. (|l|) are effective interactions mediated through several channels. The mechanisms include 
interactions between the adsorbate and the substrate electron structure, adsorbate-induced deformations of the sub- 
strate, interactions with the fluid electrolyte, and (screened) electrostatic interactions. Theoretical and computational 
methods are not yet sufficiently advanced to obtain these multi-source interactions from first principles. A method 
to estimate them on the basis of comparison between lattice-gas calculations and experiments is summarized in the 
next section. 



III. VOLTAMMETRY AND ADSORBATE STRUCTURE NEAR EQUILIBRIUM 

In the limit of vanishing potential sweep rate, the voltammetric current per unit cell is the time derivative of the 
charge transported across the interface during the adsorption/desorption process: 

q = -e{z A e A + z B B ) , (4) 

where e is the elementary charge. The current density is easily obtained in terms of the lattice-gas response functions, 
dQx/dpy, as 

P / 2 ae A do B 2 de B \dE 

A strategy has been developed to determine the effective interactions |||| , providing a practical alternative to "first- 
principles" methods. This consists in fitting thermodynamic and structural model predictions from zero-temperature 
phase diagrams and numerical simulations at room temperature to experiments, considering a variety of physical quan- 
tities. Obviously, this also has its problems. The potential number of parameters is large, and there is no guarantee 
that a minimal set of interactions is unique. Nevertheless, previous lattice-gas studies of electrochemical systems 
indicate that when all available experimental information is consistently included, this approach has considerable 
predictive power. The steps in the modeling strategy can be summarized as follows. 

1. Use prior theoretical and experimental knowledge about the substrate lattice structure and lattice constant and 
the shapes and sizes of the adsorbate particles to formulate a specific lattice-gas model. 

2. Use experimental information about adsorbate coverages and adlayer structure to determine the adsorbate 
phases. 

3. Perform a group-theoretical ground-state calculation to construct a ground-state diagram (zero-temperature 
phase diagram) and determine a minimal set of effective interactions consistent with the observed adsorbate 
phases. 

4. At nonzero temperatures, the thermodynamic and structural properties of the lattice-gas model constructed 
through steps 1-3 can be studied by various analytical and numerical methods, depending on the quantities 
of interest and the complexity of the Hamiltonian. In the applications discussed here, we have mostly used 
numerical MC simulation. This method has the advantage that it is quite accurate, even for two-dimensional 
systems ppf , and relatively straightforward to implement. 

5. The finite-temperature properties obtained in Step 4 should be used to refine the effective interactions by 
comparison with the available experiments, or by obtaining additional experimental data for such comparison. 

Steps 4 and 5 should be iterated until satisfactory agreement between model and experiment is achieved. 

Examples of the application of the approach sketched above to the chemisorption of small adsorbate particles on 
single-crystal electrodes include the adsorption of urea on Pt(100) from perchloric acid [0-|9) and the UPD of Cu on 
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Au(lll) from sulfuric acid [|L3|,[17 18|. Both systems exhibit a prominent peak sharpening in the cyclic voltammogram 
(CV) when a small concentration of the adsorbate (urea or copper ions) is added to the electrolyte. Whereas the 
urea/Pt(100) system develops a single CV peak plj , in copper UPD two peaks are exhibited [ pd|JT^ ]. The voltammetric 
changes are much weaker or absent when the same substances are adsorbed on other crystal planes of the same metals, 
indicating that they depend crucially on the geometric fit between the adsorbates and the surface. These peaks are 
associated with adsorbate phase transitions |l2|]22| l. 

Below we consider in detail simulations of some equilibrium and nonequilibrium aspects of the UPD of Cu with 
sulfate on Au(lll). 



IV. UPD OF CU WITH SULFATE ON AU(lll) 
A. Near-equilibrium aspects 

In underpotential deposition (UPD), a monolayer of one metal is electrochemically adsorbed onto another in a 
range of electrode potentials more positive than those where bulk deposition occurs. The UPD of Cu on Au(lll) 
in sulfuric acid has been intensively studied, both experimentally (see discussion of the literature in Ref. |Tt[ ) and 
theoretically [p~3| — P~S|] - The most striking feature in CV experiments on this system is the appearance of two peaks, 
100~150 mV apart, upon addition of Cu 2+ ions to the sulfuric-acid electrolyte ]ll] , fL2( |. Typical CV profiles are shown 
in Fig. |](a), together with simulation results. In the potential range between the peaks, the adsorbate layer has a 
(V3xV3) structure with 2/3 ML Cu and 1/3 ML sulfate ^3L|8|j23|, first proposed by Huckaby and Blum (HB) @. 
Typical experimental and simulated coverages are shown versus electrode potential in Fig. [l](b) . 

The lattice-gas model used by Zhang et al. [17| has interactions through fourth-nearest neighbors (see also Ref. p8|). 
Sulfate is assumed to coordinate the triangular Au(lll) surface through three of its oxygen atoms, with the fourth 
S-0 bond away from the surface. This gives a triangular "footprint" that reasonably matches the Au(lll) unit cell. 
The adsorption sites for the Cu and sulfate are assumed to lie on the same triangular lattice. The model is illustrated 
in Fig. g. 

Adsorption isotherms and CV currents at room temperature were obtained from MC heat-bath simulations on LxL 
triangular lattices with L=30 and 45. The simulated quantities were adjusted to obtain overall agreement with the 



experimental data as described in Sec. III. The experimental and simulated CV currents and coverages are shown 
together in Figs. |l|(a) and [j](b), respectively. The resulting estimates for the effective lateral interactions are given in 
the caption of Fig. g. 

It has been experimentally observed p4p^ ] that sulfate is adsorbed on top of the Cu monolayer at negative 
potentials. Zhang et al. used a simple mean-field approximation for the sulfate coverage in this second layer: 

(2) 1 

6^ = a9 c (± - 6 S ), where a is a phenomenological constant. Since the transfer of sulfate between the gold 
and copper surfaces does not involve an oxidation/reduction process, the total charge transport per unit cell during 
the adsorption/desorption process becomes 

g = - e [z s (e s + e^ ) ) + ^ce c ] , (6) 

which modifies the CV current given in Eq. (|^) in a straightforward way |T7[ . The effective electrovalences, zg and 
zq, must be determined from experiments. Zhang et al. obtained the values, zc~+l-7 and zs~— 1.1 [^7| , using data 
from Omar et al. pfij ]. To within the statistical errors, these estimates agree with independent experimental results 
by Shi and Lipkowski [^J,^5|. Comparing with values proposed by Blum et al. Jl8| based on a theoretical mean-field 
calculation, we observe that their estimate for zq is about 35% less positive than those obtained by Zhang et al. and 
Shi and Lipkowski, while there is no statistically significant disagreement about zs- 

The ground-state diagram corresponding to the interactions obtained in Ref. |l7[] is shown in Fig. ||. For large 
negative /lg, only copper adsorption is possible. Similarly, in the limit of large positive /Is and large negative the 
zero-temperature phase is the (v^x v^c/ 3 sulfate phase [we denote phases as (1x7)*] characteristic of the hard- 
hexagon model corresponding to the infinitely repulsive nearest-neighbor sulfate-sulfate interaction ^gg* ](|[l3] fl5| . 
(See Fig. |) 

The potential scan path corresponding to the CV profile and coverages shown in Fig. [j] is indicated by the dotted 
line labeled "1" in the ground-state diagram, Fig. g. Starting from the negative end, we consider a scan proceeding in 
the direction of positive electrode potential (upper left to lower right in Fig. [|). Near the CV peak at approximately 
200 mV vs. Ag/AgCl (labeled Peak # 2 in Fig. |l|), the sulfate begins to compete with copper for the adsorption 
sites, resulting in a third of the copper desorbing and being replaced by sulfate. Due to the strong effective attraction 
between the copper and sulfate adparticles, the mixed (y/3xy/3)y^ phase is formed (see Fig. ^). In this phase, which 
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extends through the entire potential region between the two CV peaks, Cu forms a honeycomb lattice with a sulfate 
molecule at the center of each cell. As the CV peak at approximately 300 mV (labeled Peak # 1 in Fig. |l|) is reached, 
most of the copper is desorbed within a narrow potential range. As it is thus deprived of the stabilizing influence of 

the coadsorbed copper, the sulfate is almost completely desorbed. 

Dynamical aspects of the phase transition at Peak # 1 are discussed in Sec. [V B below. 



B. Current transients far from equilibrium 



Whereas the preceding discussion presents a near-equilibrium theory, with voltammetric currents at slow potential 
sweep rates obtained from adsorption isotherms, analysis of fast sweeps or potential-step experiments requires a true 
nonequilibrium treatment. In the early-time regime, mean-field rate equations can be satisfactory P^^q| . However, 
for later times, when coalescence of adsorbate islands may become important, this approach is no longer necessarily 
reliable [j29|-j31[. For systems in which the phases involved are ordered, the microscopic adlayer structure and the 
dynamical details of the adsorption and lateral diffusion processes become important, even at early times. We have 
therefore initiated a study of the dynamics of electrochemical adsorption by numerical simulation of microscopic 
models. Here we present some preliminary results of this work. 



We consider the model for UPD of Cu with sulfate on Au(lll), discussed in Sec. [V A above, with the parameters 
that were obtained by Zhang et al. through comparison with near-equilibrium experiments JTt) . We focus on the 
current transients following sudden changes of the electrode potential across the discontinuous phase transition at 
Peak # 1, and we discuss our results in the light of recent experiments by Holzle et al . E9L 

The model dynamics is implemented through a refusal- free dynamic MC algorithm |3j]33j, in which the transition 
rates are determined by free energies of activation ( "free-energy barriers" ) p4[ . The transition rate from a microscopic 
configuration a to a different configuration b is taken to be 

k a ^ b = k' exp ( , (7) 



RT 

where k' is a standard rate and AG'^j is the free energy of activation for the particular transition. The eventual 
approach to equilibrium is ensured through the detailed-balance condition, 

AG^ b - AGU = H(b) - H(a) , (8) 

where 7i(a) and 7i.(b) are the values of the lattice-gas Hamiltonian in Eq. ([!]), corresponding to the two configurations. 
In these exploratory simulations we have taken the free-energy barrier associated with adsorption of particles of species 
X to depend linearly on the electrochemical potential fix (and thus on the electrode potential) through a Butler- Volmcr 
type relation @||, 

AG; dsX (/2 x ) = AG: dsX (0) - afl x , (9) 

and for simplicity we set the transfer coefficient a=l/2. The barriers associated with lateral diffusion are taken to be 
independent of the electrochemical potentials and, except insofar as they must satisfy Eq. (||), of the local geometry 
of the adsorbate layer. 

The time unit used in the simulations, called MC steps per site (MCSS), is determined by the standard rate 
constants and is expected to be proportional to the physical time. The proportionality constant must be determined 
by physico-chemical arguments and comparison with the experimentally observed timescales. 

We simulated positive-going and negative-going potential-step experiments along the scan path labeled "3" in Fig. ||, 
where Peak # 1 is centered at approximately 260 mV vs. Ag/AgCl fl7j| . Before the positive-going steps, the model 
was equilibrated in the ordered (\/3 x V%)y3 phase, 15 mV on the negative side of the transition. At time t=0 the 
potential was changed instantaneously to a final value a distance AE past the transition. Following Holzle et al., 
we used AE — 24, 28, 31, and 33 mV. For the negative-going steps we equilibrated in the disordered low-coverage 
phase, 45 mV on the positive side of the transition, followed at t—Q by a step into the ordered-phase region with 
AE = — 15 mV. Since the second-layer sulfate coverage is negligible in this potential region [see Fig. |](b)], in this 
study we only considered first-layer effects. 

The simulated current transients [in units of elementary charges per Au(lll) unit cell and MCSS] are shown in 
Fig. ||. We used systems of size 60 x 60 with periodic boundary conditions. The results were averaged over 100 trials 
for each value of AE, and the current data were smoothed using a running average over seven time points. As we 
discuss in detail below, the experimentally observed asymmetry with respect to potential steps of opposite sense is 
captured by our dynamical model. 



4 



The current transients for the positive-going potential steps, which are shown in Fig. ||(a), are characterized by 
a sharp current spike at early times, followed later by a second, broader maximum. The experimentally observed 
trends with respect to changes in AE are also reproduced, both for the locations and the sizes of the peaks. For 
the negative-going steps, no such two-peaked structure is observed. Instead, like in the experiments, a monotonically 
decreasing current is seen. This behavior is shown in Fig. |^(b). 

The surface coverages, 0c and Os, are shown versus time for positive-going steps in Figs. ||(a) and ||(b), respectively. 
Comparison of these results with Fig. |^(a) shows that the early-time current peak is almost entirely due to desorption 
of Cu. Only as 0c reaches a "metastable" value in the 0.4-0.5 ML range, does the sulfate start to desorb at a 
significant rate. At later times the relaxation appears to follow the nucleation-and-growth dynamics characteristic 
of the decay of a metastable phase jl9 27 ,3f|. The evolution of the spatial structure with time is illustrated by MC 
"snapshots" in Fig. [t]. The growth of islands of the equilibrium low-coverage phase is quite apparent. A detailed 
investigation of the time dependence in the late-time regime, including the effective "Avrami exponent" |l9|,^] and 
its dependence on the free-energy barriers, will be presented at a later date ]36| ]. 

The situation following a negative-going step is quite different. As seen in Fig. ||, Cu is rapidly adsorbed to a 
coverage of approximately 1/3 ML, which produces the large current at early times. The snapshots in Fig. || [in 
particular Fig. 0(b)] indicate that this structure is globally disordered, but locally has the symmetry of the Cu-only 
(-\/3 x Vo)i/ 3 phase. For the particular potential steps used here, the energy of this phase lies between those of the 
mixed phase and the disordered low-coverage phase. The subsequent approach towards the mixed, ordered equilibrium 
phase requires the removal of a large number of domain walls by lateral diffusion and is consequently very slow. Again, 
a detailed investigation of this dynamics and its dependence on the model parameters will be presented later |3q| . 



V. CONCLUSIONS 



In this paper we have briefly reviewed some earlier applications of statistical-mechanical lattice-gas models to 
electrochemical adsorption problems. This approach, which was originally transferred from "traditional" surface 
science, has proven to be a useful and versatile tool in the study of multicomponcnt chemisorption in electrochemical 
systems near thermodynamic equilibrium. 

Following this brief review, we presented preliminary results from a new application of the lattice-gas approach 
to dynamical phenomena far from equilibrium. In particular, we considered a dynamic extension of a model for 
UPD of Cu with sulfate on Au(lll), originally introduced by Huckaby and Blum fl3|| . We used an extension of this 
model, with parameters obtained by Zhang et al. P,|l7|] through comparison of simulation results to experiments at or 
near equilibrium. To this we added activated dynamics including adsorption, desorption, and lateral diffusion. Our 
simulated results for current transients in this system show encouraging, qualitative agreement with recent experiments 
by Holzle et al. ]19j . This agreement includes the prominent asymmetry between the transients observed in positive- 
going and negative-going potential-step experiments. Although our results leave many details for further study, we 
believe they open the door towards numerical investigations of a wide range of interesting dynamical phenomena in 
electrochemical systems. 
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FIGURE CAPTIONS 



FIG. 1. UPD of Cu with sulfate on Au(lll). Electrolyte composition: 1.0 mM CuS0 4 + 0.1 mM H2SO4. (a): CV profiles, 
AE/At = 2mV/s. Experimental (dot-dashed) and simulated (30x30: solid; 45x45: x) current densities, (b): Simulated 
coverages of Cu (dashed), first-layer sulfate (dotted), and total sulfate (solid), together with corresponding Auger Electron 
Spectroscopy data (□ and x with error bars, respectively). After Ref. JL7| . 

FIG. 2. Lattice-gas model for the UPD of copper (C) on Au(lll) in the presence of sulfate (S) p7| . The relative positions 
of copper (•) and sulfate (A) correspond to the effective interactions in Eq. ([j]). The numbers are the corresponding values of 
&$> Y , in kJ/mol. From Ref. j§. 

FIG 3. Ground-state diagram for the lattice-gas model of copper UPD on Au(lll), shown in the (fis,fic) plane. The 
effective interactions are given in Fig. [| The solid lines represent zero-temperature phase boundaries. The dotted line labeled 
"1" represents the voltammetric scan path at room temperature, corresponding to the data shown in Fig. [I]. The lines labeled 
"2" and "3" correspond to 5 mM and 0.2 mM Cu 2+ with unchanged sulfate concentration, respectively. The end points of all 
three lines correspond to _E=120 mV (upper left) and 420 mV (lower right) vs. Ag/AgCl. The solid diamonds indicate the 
positions of Peak # l ; and the solid squares indicate the positions of Peak # 2 in the simulated room-temperature CV currents. 
The phases are denoted as (AxF)| s . From Ref. 0. 

FIG 4. Ground-state configurations corresponding to the main phases in the ground-state diagram, Fig. |j. The adsorption 
sites are shown as o, and Cu and sulfate are denoted by • and A, respectively, as in Fig. |. Adapted from Ref. 0. 

FIG. 5. Simulated current transients for potential steps AE of different magnitudes and directions across the phase transition 
at Peak # 1. The potential is changed along the scan path marked "3" in Fig. and it takes the system between the ordered 
(\/3 x v^)^ phase and the disordered low-coverage phase. The time t is given in MC steps per spin (MCSS), and the current 
density in elementary charges per Au(lll) unit cell and MCSS. (a): Positive-going potential steps, AE — 24, 28, 31, and 
33 mV. (b): Negative-going potential step, AE = — 15 mV. 



FIG 6. Simulated adsorbate coverages versus time, following positive-going potential steps. The data correspond to the 
current transients shown in Fig. ||(a). (a): Copper, (b): Sulfate. 

FIG 7. "Snapshots" of configurations produced during a simulation of a 21 x 21 system following a positive-going potential 
step with AE = 24 mV. Nucleation and growth of domains of the equilibrium low-coverage phase are apparent. The times 
shown are t=0, 100, 300, and 500 MCSS. 

FIG 8. Simulated adsorbate coverages versus time, following a negative-going potential step with AE = — 15 mV. The 
data correspond to the current transient shown in Fig. @(b). Copper (solid); sulfate (dashed). 

FIG 9. "Snapshots" of configurations produced during a simulation of a 21 x 21 system following a negative-going potential 
step with AE = — 15 mV. Many locally ordered domains, separated by domain walls, are quickly formed. The subsequent 
equilibration requires the elimination of domain walls by lateral diffusion and is consequently very slow. The times shown are 
t=0, 10, 200, and 1000 MCSS. 
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